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Abstract. The emission line survey within the Calar Alto Deep Imaging Survey (CADIS) detects galaxies with 
very low continuum brightness by using an imaging Fabry-Perot interferometer. With spectroscopic follow-up 
observations of M B k, -19 CADIS galaxies using FORS2 at the VLT and DOLORES at TNG we obtained oxygen 
abundances of 5 galaxies at z ~ 0.4 and 10 galaxies at z ~ 0.64. Combining these measurements with published 
oxygen abundances of galaxies with Mb Ss — 19 we find evidence that a metallicity-luminosity relation exists at 
medium redshift, but it is displaced to lower abundances and higher luminosities compared to the metallicity- 
luminosity relation in the local universe. Comparing the observed metallicities and luminosities of galaxies at z 3 
with Pegase2 chemical evolution models we have found a favoured scenario in which the metallicity of galaxies 
increases by a factor of ~ 2 between z ~ 0.7 and today, and their luminosity decreases by ~ 0.5 — 0.9 mag. 
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1. Introduction 

In the local universe, metallicity is well correlated with the 
absolute luminosity (stellar mass) of galaxies (Skillman 
et al. 1989, Zaritsky et al. 1994, Richer & McCall 1995, 
Garnett et al. 1997, Hunter & Hoffman 1999, Melbourne 
& Salzer 2002, among others) in the sense that more lumi- 
nous galaxies tend to be more metal rich. This metallicity- 
luminosity relationship is often attributed to the action 
of galactic superwinds: massive (more luminous) galax- 
ies reach higher metallicities because they have deeper 
gravitational potentials which are better able to retain 
their gas against the building thermal pressures from su- 
pernovae, whereas low-mass systems eject large amounts 
of metal-enriched gas by supernovae-driven winds before 
high metallicities are attained (e.g., MacLow & Ferrara 
1999). However, a greater degree of gas consumption in lu- 
minous galaxies and/or suppressed infall of low-metallicity 
or pristine gas may also play a major role (Pagel 1997). 

Optical emission lines from Hn regions have long 
been the primary mean of gas-phase chemical diagnosis 
in galaxies (see, e.g., the review by Shields 1990). The ad- 
vent of large telescopes and sensitive spectrographs has 
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enabled direct measurement of the chemical properties in 
the ionized gas of cosmologically-distant galaxies with the 
same nebular analysis techniques used in local Hn regions. 
Distant galaxies subtend small angles on the sky, compa- 
rable to typical slit widths in a spectrograph. A typical 
ground-based resolution element of l'.'O corresponding to 
a linear size of ~ 6kpc at z = 0.5 encompasses the en- 
tire galaxy at medium redshift; so galaxy spectra tend 
to be integrated spectra. The use of spatially-integrated 
emission line spectroscopy for studying the chemical prop- 
erties of star-forming galaxies at earlier epochs has been 
explored by Kobulnicky et al. (1999). They concluded that 
spatially-integrated emission line spectra can reliably indi- 
cate the chemical properties of distant star-forming galax- 
ies. They also found that, given spectra with sufficient 
signal-to-noise, the oxygen abundance of these galaxies 
can be measured to within ±0.2 dex using the R23 method, 
first suggested by Pagel et al. (1979). The R23 method al- 
lows to derive metallicities also for faint galaxies where 
only prominent emission lines are observable (see more 
details in Sect. 4). 

At medium redshift, all previous attempts to study the 
evolution of metallicity with cosmic time (Kobulnicky & 
Zaritzky 1999, Hammer et al. 2001, Carollo & Lilly 2001, 
Contini et al. 2002) have been based on continuum se- 
lected samples, the emission lines of which were identified 
using spectroscopy. Owing to magnitude-limited selection 
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effects in these surveys, the samples of medium redshift 
galaxies for which metallicities were determined are bi- 
ased towards luminous galaxies. Thus, for instance, the 
CFRS sample used by Carollo & Lilly (2001), selected by 
/ < 22.2, contains only bright galaxies (M B < —20) at 
redshift 0.6 < z < 1, the continuum of which is domi- 
nated by an evolved stellar population. 

Kobulnicky & Koo (2000) found that star-forming 
Ly man-break galaxies at z ~ 3 are 2— 4 mag more lumi- 
nous than local spiral galaxies of similar metallicity, and 
thus are offset from the local luminosity-metallicity rela- 
tion. Also Pettini et al. (2001) found that Lyman break 
galaxies at z <~ 3 are significantly overluminous for their 
metallicities. 

It seems therefore possible that the whole metallicity- 
luminosity relation is displaced to lower abundances at 
high or medium redshifts, ilustrating how galaxies partic- 
ipate in the chemical evolution process. 

The Calar Alto Deep Imaging Survey (CADIS, 
Meisenheimer et al. 1998, 2004) allows the selection of 
galaxies by their emission line fluxes detected via narrow- 
band (Fabry-Perot) imaging, regardless of their contin- 
uum brightness. Thus, CADIS can reach galaxies with 
very low continuum brightness (and absolute luminosity), 

1. e. galaxies which have such faint continuum that they are 
only detectable by their emission lines. Combining metal- 
licities for such faint CADIS galaxies in several narrow red- 
shift bins with published abundances for more luminous 
galaxies at medium z we aim to establish the metallicity- 
luminosity relation for several look-back time bins in the 
range 0.3 < z < 1. The metallicity-luminosity relation at 
different reshifts is a valuable probe of chemical galaxy 
evolution, and a consistency check for chemical evolution 
models. 

This paper is structured in the following way: In Sect. 

2 we describe how we select galaxies at 0.39 z Si 0.65 
from the CADIS emission line sample in order to deter- 
mine their metallicities by follow-up spectroscopy. In Sect. 

3 we present the spectroscopic follow-up of such emission 
line galaxies. In Sect. 4 we compute oxygen abundances 
of the spectroscopic sample. Finally, in Sect. 5 we dis- 
cuss the evolution of the metallicity-luminosity relation 
with redshift. The "concordance" cosmology with H = 70 
km s _1 Mpc -1 , flo = 0.3, Q\ = 0.7 is used throughout 
this paper. For the solar oxygen abundance we use the 
value 12 + log(0/H) = 8.87 from Grevesse et al. (1996). 
Note that "metallicity" and "abundance" normally de- 
notes "oxygen abundance" throughout this paper. 

2. Selection of emission line galaxies at medium 
redshift from CADIS 

CADIS combines a moderately deep multi-band survey 
(10 a limit, Ru m = 24) with a deep emission line survey 
employing an imaging Fabry- Perot-Interferometer (Fii m = 
3 • 10~ 20 Wm~ 2 ), which covers three waveband windows 
(henceforth called FP windows) essentially free of atmo- 
spheric OH emission lines: A (A = 702 ± 6nm), B (A = 



819 ± 5 nm), and C (A = 920 ± 8 nm). An emission line de- 
tected in the Fabry-Perot scan can be any of the prominent 
nebular lines, e.g., Ly-a, [Oil] A 3727, H/3, [O III] A 5007, 
or Ha. The line identification is done by means of the 
line profile, by fitting template galaxy spectra to the ob- 
served multi- filter SED, and by considering the flux in 
"veto-filters" at the expected wavelengths of other lines 
(Hippelcin et al. 2003, Meisenheimer et al 2004). 

The analysis of FP window A and B is nearly com- 
pleted in four of the 6 CADIS fields with FP observa- 
tions (Table 1 in Maier et al. 2003). In these four CADIS 
fields we have found 128 galaxies by their [O III] A 5007 
emission detected in the FP windows A or B; they have 
redshifts in the redshift bin 0.392 < z < 0.415, and 
0.625 < z < 0.648, respectively. Figures 1 and 2 show 
two examples of galaxies detected by their [O III] A 5007 
emission line in the FP scan. The [O III] A 5007 line of the 
galaxy 23h-671455 (Fig. 1) is detected in the FP window 
A, and we see the [O II] A 3727 emission line in the veto- 
filter at 522 nm. Moreover, Ha and the [N II] AA 6548, 6584 
doublet are seen in window C. For the galaxy Olh-584247 
(Fig. 2), the [O III] A 5007 line is seen by the Fabry-Perot 
in window B, and we detect the [Oil] A 3727 line in the 
veto-filter at 611 nm. 

From the sample of 128 objects detected by their 
[O III] A 5007 line in FP windows A or B we selected 
emission line galaxies for spectroscopic-follow up in 
order to determine their metal abundances. Oxygen 
lines are the most useful for measuring the abun- 
dance, since all important ionization stages can be ob- 
served. The temperature-sensitive [O III] A 4363 line is of- 
ten too weak to be measured even in the local universe 
([OIII] A 4363/ [O III] A 5007 ^ 1/50). To overcome this 
problem, Pagel et al. (1979) suggested the empirical abun- 
dance indicator 

I([OII] A3727) +I([OIII] AA4959,5007) 

R23 im { ) 

for faint objects, such as galaxies at high redshift. 
The method has been later refined by McGaugh (1991) 
using a set of photoionisation models, and observation- 
ally confirmed to be reliable by Kobulnicky et al. (1999). 
For the 128 selected objects the [OIII] A 5007 line is 
measured in the respective Fabry-Perot window, while 
H/3 must be measured by spectroscopic follow-up, and 
[Oil] A 3727 is detected by a veto filter at 522 or 611 nm. 
However, since the determination of the [O II] A 3727 flux 
relies only on a medium band filter, this flux is often 
not accurate enough to determine oxygen abundances 
with the R23 relation. Therefore, spectroscopic follow-up 
is necessary to measure not only the H/3 line, but also 
[Oil] A 3727. Since the [O III] AA 4959, 5007 and H/3 lines 
lie close in wavelength to each other, the accurate fluxes 
of the [OIII] AA 4959, 5007 lines are provided for free by 
the spectroscopic follow-up. 

From our 128 objects detected by their [OIII] A 5007 
emission line seen in the Fabry-Perot we obtained spec- 
tra of 21 objects with absolute magnitudes fainter than 
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Fig. 1. The galaxy 23h-671455 at z » 0.41: Photometry in all 14 optical CADIS filters fitted by a continuum-model with 
overlayed prominent emission lines (left panel), the Fabry-Perot measurements in window A with a [OIII] doublet profile fitted 
to the observed flux data (center), and the Fabry-Perot measurements in window C with a Ha + [Nil] AA6548, 6584 three-line 
profile fitted to the observed flux data (right panel). Note that the [Oil] A 3727 line is seen in the veto- filter at 522 nm. The 
numbers in the lower right edges of the right panels specify the signal-to-noise ratios, Fi ine / ai ine ; Fu ne and ou ne are derived by 
fitting the Fabry-Perot instrument profile to the observed flux data points; see also Hippelein et al. (2003). 




Wavelength [nm] 

Fig. 2. The galaxy Olh-584247 at z » 0.64: Photometry in all 14 optical CADIS filters fitted by a continuum-model with 
overlayed prominent emission lines (left panel), and the Fabry-Perot measurements in window B with a [OIII] doublet profile 
fitted to the observed flux data (right panel). Note that the [Oil] A 3727 line is seen in the veto-filter at 611 nm. 



Mb ~ —19, which were able to be placed on the slitmasks 
together with Ly-a candidates (Maier et al. 2003) in order 
to maximize the number of galaxies observed; see Section 
3 for the details of the observations. The large wavelength 
coverage of the CADIS filters allows to derive precise ab- 
solute magnitudes without using any K-correction factor, 
since the flux value at the restframe wavelength is calcu- 
lated by linear interpolation between the measurements 
of two adjacent filters. In order to extend the sample 
of already published metallicities of galaxies at medium 
redshift to fainter absolute magnitudes, we selected for 
the spectroscopic follow-up objects with Mb ^ —19. We 
removed five objects from this sample because the H/3 
line had too low signal-to-noise ratio for reliable oxygen 
abundance determinations (see Kobulnicky et al. 1999 
for the discussion of the signal-to-noise required in or- 
der to get reliable metallicity measurements with the R23 
method). The remaining 16 galaxies are listed in Table 



2. As shown in Fig. 3, this final sample is only slightly 
biased towards special F([OIH] A 5007) /F( [Oil] A3727) 
ratios a ratio which decreases with metallicity 
(Stasinska & Leitherer 1996) - in the sense that objects 
with F([OIII] A5007)/F([OII] A3727) < 1 are underrepre- 
sented (28% in the full sample, 12% in the spectroscopic 
sample). 

3. Spectroscopic observations 

The spectroscopic follow-up observations of emission line 
galaxies in two redshift bins, at 0.392 ^ z ^ 0.415, and 
at 0.625 < z < 0.648, in the CADIS Olh- and 23h- 
fields were obtained in the summer and autumn of 2001, 
and in the summer of 2002, using FORS2 at the VLT. 
Galaxies in 09h-field, and one galaxy in Olh-field at z as 
0.64 (Olh-5085) were observed with the Low Resolution 
Spectrograph DOLORES at Telescopio Nazionale Galileo 
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Fig. 3. The histogram of F([0 III] A 5007)/F[O II] A 3727) flux 
ratios (as determined by the CADIS filter and Fabry-Perot 
measurements) of 128 emission line galaxies selected by their 
[O III] A 5007 line seen in the Fabry-Perot; the filled boxes show 
the CADIS flux ratios for the galaxies with follow-up spec- 
troscopy listed in Table 2. Note that the negative values (dotted 
lines) come from the fact that, for galaxies with undetectable 
[O II] A 3727 flux, the uncertainties in the interpolated contin- 
uum and the background noise place the formally measured 
F([OII] A 3727) flux symmetrically around F([Q II] A 3727)= 0. 



(TNG). The slitmasks contained l'.'O up to l'/6 wide slits; 
the length of the slits varied between lO'.'O, and 2f/.'0. Since 
we found from our first observations that the l'.'O slit may 
miss (part of) the line emitting region (see discussion in 
Maier et al. 2003), we used wider slits (up to 1'.'6) for the 
observations with TNG (masks e, /), and for the VLT 
observations in summer 2002 (masks d, i). 

For the VLT observations we used three different 
grisms, the lower resolution 3001 grism, which gives a 
spectral resolution of about 1.2 nm at 800 nm for l'.'O wide 
slitlets; the 600 RI grism, which gives a higher spectral res- 
olution of about 0.8 nm at 800 nm for l'.'O wide slitlets, and 
the 600R grism, which gives a spectral resolution of about 
0.5 nm FWHM at 600 nm. However, since the objects have 
to be put on the masks according to the wavelength range 
to be observed, the 600RI grism has a smaller spectro- 
scopically field-of-view than the 3001 grism, and thus al- 
lows less freedom in target positions. Moreover, one can 
measure only the [O III] AA 4959, 5007 and H/3 lines of a 
galaxy at z ~ 0.64, and, in order to get [O II] A 3727, the 
600R grism has to be used additionally. 

Spectra were reduced using the package LONG pro- 
vided by MIDAS. Images were bias-subtracted first. In 
order to remove the low spatial frequencies along the dis- 
persion axis from the flatfield, we averaged the original 
image along the slit, fitted the resulting one-dimensional 
image by a polynomial of second order, and divided the 
original flatfield by this image to obtain a normalized flat- 
field. Every mask contained about 20 single spectra. Every 
single spectra was extracted, treated as a long slit spectra, 



and divided through the respective normalized flatfield. 
For the wavelength calibration night sky lines were used, 
which have the advantage to be taken with the telescope 
and spectrograph in the same orientation as the science 
spectra. The one-dimensional spectra of each object were 
extracted with an aperture of about 10 pixels, i.e., about 
2", using the algorithm by Home (1986). Fluxes were cal- 
ibrated from digital units (ADUs) to physical flux units 
(W nr 2 s- 1 A" 1 ) using multiple observations of the spec- 
trophotometric standard stars LTT 7379, EG 274, and 
LTT 7987 (Hamuy et al. 1992, 1994). 

For the observations with TNG we used the HRR and 
MRB grisms for galaxies at z « 0.4, and the LRR grism 
for the galaxy Olh-5085 at z ~ 0.64. Data reduction of 
the spectra was performed similar to the VLT data re- 
duction. For the flux calibration the spectrophotometric 
standard star HD93521 (Oke 1990) was used. Table 1 lists 
the details of the observations with the VLT and TNG. 

4. Flux ratios and oxygen abundances 

4.1. Emission line measurements 

Emission line fluxes were measured interactively using the 
integrate/line routine in MIDAS. The flux error is calcu- 
lated as x (j, where N is the number of pixels summed 
in a given emission line, and a is the root-mean-squared 
variation in a region of the spectrum next to the emis- 
sion line. This yields typical errors of 5% for [O III] A 5007, 
[O III] A 4959, and [Oil] A 3727, and 5% to 12 % for H/3, 
depending on the position of H/3 between the night sky 
lines.. The error of R23 is thus dominated by H/3. Flux ra- 
tios were corrected for underlying Balmer absorption with 
a general H/3 absorption equivalent width of 2A , represen- 
tative of local irregular, HII, and spiral galaxies (see e.g., 
McCall et al. 1985, Skillman & Kcnnicutt 1993, Izotov et 
al. 1994), and in agreement with the newest results for R(3 
absorption equivalent widths from the study of about 600 
SDSS strong emisison line galaxies (Kniazev et al. 2004). 

4.2. Extinction 

The dereddened value for a flux line ratio, I(A1)/I(A2), is 
given by: 



I(A1) _ F(A1) 
I(A2) ~ F(A2) 



x 10 c(/(Al)-/(A2)) 



(2) 



where F(A) is the observed flux at a given wavelength, 
c is the logarithmic reddening parameter that describes 
the amount of reddening relative to H/3, and /(A) is 
the wavelength-dependent reddening function (Whitford 
1958). The value of c can be extimated from the relation 
c = 1.47 x E B -v from Scaton (1979). 

The ratio I([0 III] A 5007)/I(H/3) which goes into the 
R23 relation is: 



I([OIII] A 5007) _ F([OIII] A 5007) 



I(H/3) 



F(H/J) 



x 10 



-0.034 c 



(3) 
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Field 


Mask 


t(sec) 


Grism 


Res.(nm/slit) 


Slitwidth 


Inst. /Telescope 


Olh 


a 


12860 


3001 


1.2 


1".0 


FORS2/VLT 


Olh 


b 


12500 


600RI 


0.8 


1".0 


FORS2/VLT 


Olh 


b 


7500 


600R 


0.5 


1".0 


FORS2/VLT 


Olh 


c 


6085 


600RI 


0.8 


1".0 


FORS2/VLT 


Olh 


c 


1700 


600R 


0.5 


1".0 


FORS2/VLT 


Olh 


d 


15750 


3001 


1.7 


1".4 


FORS2/VLT 


Olh 


e 


6000 


LRR 


1.2 


l".l 


DOLORES/TNG 


09h 


f 


15600 


HRR 


0.5 


l".6 


DOLORES/TNG 


09h 


f 


3000 


MRB 


1.0 


1".6 


DOLORES/TNG 


23h 


g 


16080 


3001 


1.2 


1".0 


FORS2/VLT 


23h 


h 


9000 


600RI 


1.8 


1".0 


FORS2/VLT 


23h 


i 


15750 


3001 


1.7 


1".4 


FORS2/VLT 



Thus, since the [O III] AA 4959, 5007 and H/3 lines lie 
close in wavelength to each other, the extinction is small 
and can be neglected for the I([0 III] A 5007)/I(H/3) ratio 
(the same applies also to the I([0 III] A 4959)/I(H/3) ratio). 

The other ratio, I([0 II] A 3727)/I(H/3), which goes into 
the R23 relation is 



I([OII] A 3727) _ F([OII] A 3727) 



W) 



F(H/3) 
F([0 II] A 3727) 



x 10 



0.308 c 



x a, 



(4) 



with a = io°- 308c . Thus, the effect of reddening in the 
R23 expression can be written as: 



R 



23 



a x F([0 II] A 3727) + F([0 III] AA 4959, 5007) 



(5) 



Assuming case B Balmer recombination, with a 
temperature of 10 000 K, and a density of 100 cm~ 3 
(Brocklehust 1971), the predicted dereddened intensity 
ratio of Ha to H/3 is 2.86 (Ostcrbrock 1989). The ef- 
fect of reddening on the ratio Ha/H/3 can be written 
as I(Ha)/I(H/3) = F(Ha)/F(H/3) x lo(-°- 332 ) c . Using the 
spectroscopic follow-up measurements of the Ha line for 
three galaxies atzw 0.4 (see, e.g., Fig. 5), and the CADIS 
Fabry-Perot measurements of the Ha line available in win- 
dow C for galaxies atz« 0.4 (see Fig. 1), the values of 
c for the 5 galaxies at z ~ 0.4 are of order of — 0.1, 
yielding a value of a in the range 1 to 1.07. This means 
that the extinction (see equation 5) can be neglected with 
respect to the errors for the faint emission line galaxies 
(Mb ^ —19) at z m 0.4, when calculating R23. We assume 
that the same is true for galaxies with Mb ^ —19 in the 
redshift bin 0.625 < z < 0.645. 

This assumption is supported by the fact that the 
[O III] A 5007 line flux is larger than the [O II] A 3727 line 
flux for the 11 galaxies atzw 0.64 of our sample (see Table 
2). Therefore, even allowing for a wider range 1 < a < 2 
and typical large F([0 III] A 5007)/F([O II] A 3727) ratios, 
equation 5 shows that the [O III] AA 4959, 5007 and H/3 



measurement errors, and not the unknown value of a, 
dominate the error of the R23 ratio. Therefore, we do not 
correct the line fluxes for extinction for the galaxies in our 
sample. 

4.3. The R23 relation 

Because the temperature-sensitive line [O III] A 4363 is not 
detectable in the spectra of galaxies at z ^ 0.4, the oxy- 
gen abundance has to be determined using the empiri- 
cal calibration between oxygen abundance and R23. Pagel 
et al. (1979) noted, based on a sample of extragalactic 
HII regions, that the measured T e , O/H and R23 were 
all correlated. This works because of the relationship be- 
tween O/H and nebular cooling: For metal-rich regions, 
the cooling in the ionized gas is dominated by emission 
in IR fine-structure lines (primarily the [O III] 52 /im and 
88 /im lines), so as O/H increases, the nebula becomes 
cooler. In response to that, the highly excited optical for- 
bidden lines, the [OIII] lines, become weaker as O/H in- 
creases (excitation goes down as T decreases). At very 
high abundances, 12 + log(0/H) ^ 8.7, H11 regions are 
very cold (T e < 6000 K) and values of R23 are small. For 
12+log (O/H) <, 8.2, the relation between R 23 and O/H 
reverses, such that R23 decreases with decreasing abun- 
dance. This occurs because at very low metallicities the 
IR fine-structure lines no longer dominate the cooling be- 
cause of the lack of heavy elements. As a result, the for- 
bidden lines reflect the abundance in the gas almost in a 
proportional way. 

As a consequence, the R23 indicator is not a mono- 
tonic function of oxygen abundance (see Figure 4). At 
a given value of R23, there are two possible choices of 
the oxygen abundance. Roughly, the low abundance or 
"lower branch" is defined by 12 + log(0/H) ^ 8.2, whereas 
the high abundance or "upper branch" is defined by 
12 + log(0/H) > 8.5. Therefore, an additional indica- 
tor (e.g., the [Nil] A 6584 line) is needed in order to re- 
solve the degeneracy in R23. About 50% of the emis- 
sion line galaxies in our sample have logR23 m the range 
0.8 — 1 placing these galaxies in the "turnaround" region 
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plication of the R23 method. The division between the 
upper and the lower branch of the R23 relation oc- 
curs around [Nil] A 6584 /Ha - 0.1. [Nil] A6584/Ha 
line ratios were measured by follow-up spectroscopy for 
three CADIS galaxies at z « 0.4. 09h-448495 has an 
[N II] A 6584 /Ha line ratio > 0.1, placing this galaxy 
on the upper branch of the R23 relation; the galaxy 
23h-671455 has a [Nil] A6584/Ha line ratio of - 0.1 
and lies in the turnaround region; and 23h-683506 has 
a [Nil] A 6584 /Ha line ratio < 0.1, placing this galaxy 
on the lower branch of the R23 relation (see spectrum in 
Figure 5, lower panel). 



Fig. 4. The oxygen abundance 12+log(0/H) as a function < 
the line ratio logR23- The calibration between R23 and O/] 
(solid lines) using the models from McGaugh (1991) shows th 
effect of varying the ionization parameter in terms of the ot 
servable line ratio [O III] A A 4959, 5007/ [Oil] A 3727 . On th 
lower, metal-poor branch, the ionization parameter become 
important. The ratio [O III] /[Oil] of 10, 1.0, and 0.1 corn 
spond (very roughly) to ionization parameters U of 10 _1 , 10 - ' 
and 10~ 4 . The two horizontal dotted lines separate the diffei 
ent regions of oxygen abundances. The arrow and the hoi 
izontal dashed line indicate the solar oxygen abundance 
12 + log(0/H) = 8.87 (Grevesse et al. 1996). 



0) 



x 
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8.2 < 12 + log(0/H) < 8.5 of the oxygen abundance vei 
sus R23 diagram (see Fig. 4). For the remaining galaxk 
we describe in the following how the decision, whether 
galaxy lies on the upper or on the lower branch, has bee" 
taken. 



4.3.1. Breaking the R23 degeneracy by the 
[Nil] A6584/Ha ratio 

The [Nil] A6584/Ha line ratio can be used to brea 
the degeneracy of the R23 relation, if one is able t 
measure and separate [Nil] A 6584 from Ha. Denicolo 6 
al. (2002) calibrated the N2 estimator, defined as N2 = 
log(I([NII] A6584)/I(Ha)), vs. the oxygen abundance, us 
ing a sample of Hn galaxies having accurate oxygen abur 
dances, plus photoionization models covering a wide rang 
of abundances. 

When the secondary production of nitrogen don 
inates, at somewhat higher mctallicity, the line rati 
[Nil] A6584/Ha increases with oxygen abundance. Au 
very low mctallicity, N2 scales simply as the nitrogen 
abundance to first order. However, in this mctallicity 
regime, the nitrogen abundance shows a large scatter rel- 
ative to the oxygen abundance, since the nitrogen abun- 
dance is much more sensitive to the history of star for- 
mation in the galaxy considered. As a result, N2 is prob- 
ably not very useful to estimate oxygen abundance ex- 
cept as a means of determining the branch for the ap- 
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Fig. 6. VLT Spectra using the 3001 grism with FORS2 of two 
galaxies at 2 « 0.64: 23h-5159 with Mb = —17.5 (upper panel), 
and 23h-6512 with M B = -18.3 (lower panel). The "features" 
observed in the continuum are largely due to residuals from 
the subtraction of strong OH emission lines. 
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Fig. 5. VLT Spectra (mask i) using the 3001 grism with FORS2. The "features" observed in the continuum are largely due to 
residuals from the subtraction of strong OH emission lines. 

Upper panel: The galaxy 23h-671455 (Mb = -19.0) at z w 0.4 shows a [Nil] A 6584 /Ha line ratio of ~ 0.1, and lies in the 
turnaround region of the R23 relation. 

Lower panel: The galaxy 23h-683506 (Mb = -17.3) at z ^ 0.4 shows a [Nil] A 6584 /Ha line ratio < 0.1, placing this galaxy 
on the lower branch of the R23 relation. 



4.3.2. Breaking the R23 degeneracy by high 
[O III] A 5007/ [Oil] A 3727 flux ratios 

From the remaining galaxies without measured Ha and 
[N II] A 6584 (which would require near- infrared spec- 
troscopy for z — 0.64) three galaxies, 09h-542442, Olh- 
246610, and 23h-408246, show a high [O III] A 5007 to 
[O II] A 3727 flux ratio (greater than five). If we put these 
galaxies on the upper branch, we would get metallicities 
of 0.5 Zq < Z < IZq. There are no galaxies in the lo- 
cal universe in this high mctallicity range which show a 
such high [O III] A 5007/ [O II] A 3727 ratio. Assuming that 
the physical properties of the interstellar medium are the 
same in the local universe and at medium redshift, these 
galaxies cannot have a such high metallicity. Therefore, 
we have to put them on the lower branch. We used this 
criterion also for the galaxies 23h-441445 and 23h-552534 
(spectra shown in Figure 6). 

We exclude the galaxy 01h-628521 oxygen abun- 
dance (indicated in Table 2 as (?)) from the follow- 
ing discussion of oxygen abundances, since we cannot 
be sure on which branch of the R23 relation to place 
it ([OUT] A 5007/ [Oil] A 3727 ~ 1). Near infrared spec- 
troscopy will be required in order to measure the Ha 
and [N II] A 6584 lines of this galaxy, and to determine 
on which branch of the R23 relation it has to be placed. 

4.4. Oxygen abundances 

We use the most recent R23 formulation by Kobulnicky et 
al. (1999) of the analytical expressions by McGaugh (1991) 



to determine the oxygen abundances of the emission line 
galaxies in our sample. These formulae express the oxygen 
abundance, 12+log (O/H), in terms of R23 and the ioniza- 
tion index [O III] AA 4959, 5007/ [Oil] A 3727. The branch 
of the R23 relation was selected as described in section 4.3. 
The computed oxygen abundances are shown in Table 2. 
Additional to the measurement uncertainties of the emis- 
sion line ratios, an uncertainty of <~ 0.12 dcx in the oxy- 
gen abundance has been assumed, which comes from un- 
certainties in the photoionization models and ionization 
parameter corrections for the R23 method. 

5. The metallicity-luminosity relation at medium 
redshift 

Oxygen abundances have been determined for a total 
number of 15 CADIS emission line galaxies at medium 
redshift with faint absolute magnitudes (Mb > —19). 
Combining these abundances with published results for 
galaxies with brighter absolute magnitudes we can study 
the metallicity-luminosity relation over a large range of 
luminosities at look-back times of 4.3Gyrs (z w 0.4), and 
5.9Gyrs (z w 0.64). Fig. 7 shows the oxygen abundance 
as a function of absolute magnitude for two redshift bins, 
at z w 0.4, and z w 0.64, respectively. Additional to the 
CADIS measurements, oxygen abundances from literature 
for galaxies in the two redshift bins are shown. It is ob- 
vious from the diagrams that only the CADIS emission 
line sample provides a meaningful sample of Mb > —19 
galaxies to study the faint absolute magnitudes part of the 
metallicity-luminosity relation at medium redshift. 
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Table 2. Emission line galaxies with unreddened line strengths and computed R23 ratios. a to 1 indicates the mask with which 
the object was observed (see Table 1). Column Branch indicates if the object lies on the lower (L) or on the upper (U) branch of 
the O/H versus R23 relationship, falls in the turnaround region (T), or if no decision between the upper (U) and lower branch 
can be made (?). The table is ordered in bins of increasing redshift, at 0.392 < 2 < 0.415, at 0.625 < z < 0.648, and z > 0.68. 



# 


Zspec 


Mb 


[OIII]/[OII] 


[OIII]/H/3 


logR,23 


Branch 


0/H L 


O/Hu 


12+log(0/H) 


09h-542442- / 


0.3985 


-18.1 


11.11 ±0.79 


3.03 ± 0.34 


0.77 ±0.10 


L 


7.64 


8.74 


7.64 ±0.15 


09h-448495 / 


0.4080 


-19.2 


1.49 ±0.11 


2.38 ±0.27 


0.69 ±0.10 


U 


7.81 


8.72 


8.72 ±0.15 


09h-392705 / 


0.4090 


-18.4 


0.91 ±0.06 


2.94 ±0.33 


0.85 ±0.10 


T 


8.16 


8.51 


8.34 ±0.15 


23h-671455* 


0.4068 


-19.0 


0.41 ±0.07 


1.96 ±0.13 


0.87 ±0.05 


T 


8.30 


8.44 


8.37 ±0.15 


2311-683506* 


0.4065 


-17.3 


1.95 ±0.61 


3.05 ±0.23 


0.77 ±0.05 


L 


7.90 


8.65 


7.90 ±0.15 


01h-628521 d 


0.6370 


-18.7 


1.23 ±0.07 


2.30 ±0.22 


0.71 ±0.02 


7 


7.86 


8.70 


7 


01h-584247 c 


0.6358 


-19.2 


1.59 ±0.11 


4.03 ± 0.45 


0.89 ±0.10 


T 


8.11 


8.51 


8.31 ±0.15 


01h-537252 c 


0.6372 


-19.1 


1.25 ±0.09 


3.70 ±0.41 


0.90 ±0.10 


T 


8.16 


8.48 


8.32 ±0.16 


Olh-471561 a 


0.6352 


-18.6 


2.33 ±0.17 


3.70 ±0.41 


0.84 ±0.09 


T 


7.97 


8.60 


8.29 ±0.31 


01h-246610 i ' 


0.6390 


-17.3 


12.50 ± 1.45 


4.76 ±0.53 


0.83 ±0.10 


L 


7.77 


8.67 


7.77 ±0.15 


01h-177719 e 


0.6221 


-18.4 


2.86 ±0.20 


5.03 ±0.56 


0.95 ±0.08 


T 


8.19 


8.46 


8.33 ±0.14 


2311-441445* 


0.6435 


-17.5 


1.91 ±0.21 


1.67 ±0.16 


0.53 ±0.03 


L 


7.53 


8.85 


7.53 ±0.15 


23h-558487* 


0.6433 


-18.4 


2.22 ±0.27 


4.00 ± 0.45 


0.86 ±0.02 


T 


8.05 


8.55 


8.30 ±0.15 


23h-552534* 


0.6447 


-18.3 


2.46 ±0.25 


2.16 ±0.20 


0.59 ±0.03 


L 


7.58 


8.82 


7.58 ±0.15 


23h-408246 9 


0.6390 


-18.6 


6.67 ±0.47 


5.04 ±0.56 


0.85 ±0.09 


L 


7.84 


8.65 


7.84 ±0.15 


23h-431392 9 


0.6440 


-19.2 


1.11 ±0.08 


3.33 ±0.37 


0.87 ±0.08 


T 


8.12 


8.57 


8.35 ±0.23 


23h-385589 ! 


0.6880 


-18.9 


0.50 ± 0.06 


2.25 ±0.25 


0.88 ±0.03 


T 


8.29 


8.44 


8.32 ±0.30 


23h-645522* 


0.8903 


-19.5 


0.87 ±0.08 


2.34 ±0.23 


0.77 ±0.02 


T 


8.03 


8.61 


8.37 ±0.15 



Since our number statistics of oxygen abundances at 
medium redshift are limited, we applied a linear least 
squares fit to the data. There may be some curvature 
in this correlation. However, it is not possible to testify 
such a curvature with our limited sample. For emission 
line galaxies at z « 0.4 we find a metallicity-luminosity 
relation: 

12 ± log(0/H) = (4.59 ± 0.91) - (0.199 ± 0.046)M B , (6) 

with an RMS scatter of 0.219. For emission line galaxies 
at z ~ 0.64 we find a metallicity-luminosity relation: 

12 ± log(0/H) = (2.81 ± 0.81) - (0.288 ± 0.042)M B , (7) 

with an RMS scatter of 0.225. The dotted (z ~ 0.4) and 
dashed-dotted line (z ~ 0.64) in Figure 7 are these linear 
fits showing that a metallicity-luminosity relation exists 
at medium redshift over a larger luminosity range. 

Fig. 8 compares the observed metallicity-luminosity re- 
lations at z = 0.4 and z = 0.64 (derived in Fig. 7) with 
the metallicity-luminosity relation for z — emission line 
galaxies (Melbourne & Salzer 2002) and the oxygen abun- 
dances of luminous galaxies at z <~ 3.1 (Kobulnicky & Koo 
2000, Pettini et al. 2001). Moreover, we compare the ob- 
served mctallicities and luminosities with Pegase2 (Fioc 
& Rocca-Volmerange 1999) models. 

5.1. Comparison to Pegase2 Models 

Pegase2 is a galaxy evolution code which allows to study 
galaxies by evolutionary synthesis. Here is a discussion of 
the parameters which were chosen by us: 



— ip(IMF): The shape of the stellar initial mass function. 
We adopt the Salpeter value, a = —2.35, between 0.1 
and 120 M (like Baldry et al. 2002, e.g.). 

— Y: The chemical yields from nucleosynthesis. We use 
Woosley & Weaver 1995 B-series models for massive 
stars. 

— M to t- The total mass of gas available to form the 
galaxy. We chose three representative masses of 2 x 
10 10 Mq, 6 x 10 10 M Q , and 2 x 10 11 M Q . 

— ti n f a u: The timescale on which the galaxy is assembled. 
We assume that galaxies are built by continuous infall 
of primordial gas (zero metallicity) with an infall rate 
that declines exponentially, ti n f a u, as implemented in 
Pegase2, 

P -t/t infaU 

M(t) = M tot — . (8) 

tin fall 

We chose two representative gas infall timescales of 
tin fail = 1 Gyr and t lnfaU = 4 Gyrs. 

— xjj(SFR): The form of the star formation rate. We 
adopt an exponentially decreasing star formation rate, 
as implemented in Pegase2: 

e -t/ P i 

^(t)= p 2— -. (9) 
pi 

— Ay: Extinction due to dust. An inclination- averaged 
extinction prescription as implemented in Pegase2 was 
included, but extinction does not play an important 
role in the present comparison (it changes the model 
B magnitudes by only 0.2 mag). 

We assume that galaxies begin to form about 1 Gyr 
after the Big Bang so that local z = galaxies have 
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Fig. 7. The oxygen abundance as a function of absolute magnitude, for z ~ 0.4 (upper panel) and z ~ 0.64 (lower panel). 
The CADIS metallicity measurements are shown as filled squares, open squares show the abundances found by Kobulnicky 
& Zaritsky (1999) at z « 0.4, and by Carollo & Lilly (2001) at z « 0.64, for more luminous galaxies. Open triangles denote 
measurements from Kobulnicky et al. (2003). The dotted (z ~ 0.4) and dashed-dotted line (z ~ 0.64) are linear fits to the data 
points. 



ages of about 13 Gyrs (in a 13.5 Gyr old universe). 
We explored a range of SFR, varying pi and p2, and 
the timescale on which the galaxy is assembled, t, n / M- 
Assuming pi — t in f a u a high metallicity is reached after 
less than 1 Gyr, and the metallicity do not change much 
in the next 12 Gyrs, not reproducing the observed abun- 
dances. Models with pi < ti n j a u are unphysical, since 
SF ceases, although gas is still infalling. Therefore we 
explored different values of pi > Unfall- It turned out 
that, in order to explain the metallicity- luminosity evo- 
lution of both luminous galaxies at z ~ 3 and galaxies 
at medium redshift and z = 0, (representative) mod- 
els with values of (t in f aU = 4 Gyrs, pi = 8 Gyrs) or 



{Unfall = 1 Gyrs, pi — 8 Gyrs), respectively, can repro- 
duce the observed metallicities. The values of p2 has been 
set to p2 = 0.2M Q ,p2 = 0.5M Q and p2 = 1M Q according 
to the total mass of a galaxy of 2 x 10 10 M Q , 6 x 10 10 M , 
and 2 x 10 11 M , respectively. Thus, only a SFR history 
with a long exponential timescale (~ 8 Gyrs) can repro- 
duce the observations, but the gas infall has to cease after 
1-4 Gyrs. 

The internal scatter of the metallicity-luminosity rela- 
tion is relatively well understood in the context of a burst- 
ing star formation mode, shown in Mouhcine & Contini 
(2002). Therefore, we applied also such a star formation 
history in our Pegase2 models: we investigated bursting 
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Fig. 8. The metallicity-luminosity relation over a look-back time of ~ 12 Gyrs. The solid line shows the luminosity-metallicity 
relation found by Melbourne & Salzer (2002) for a large sample of local emission-line galaxies (look back time OGyrs). Also 
shown are the metallicity-luminosity relations for galaxies in at z ~ 0.4 (look back time 4.3 Gyrs), and z ~ 0.64 (look back 
time 5.9 Gyrs) from Fig. 7. The location of bright galaxies at z ~ 3.1 is shown as a box encompassing the range of O/H and 
Mb derived for these objects (from Fig. 7 in Pettini et al. (2001), which takes also into account the metallicities derived by 
Kobulnicky & Koo (2000) for z ~ 3.1 objects). The dashed curves and symbols track the chemical and luminous evolution of 
PEGASE2 model galaxies with an exponential star formation rate. The galaxies are build by exponentially-decreasing infall of 
primordial gas (zero metallicity) with infall timescales of lGyr and 4Gyrs, respectively. The symbols show galaxies at a look 
back time of 0, 1, 5, 9, 11 and 12 Gyrs, assuming that the galaxies begin to assemble at a look back time of 13Gyrs. 



star formation models assuming several bursts of 50 Myrs 
duration with an inter-burst period of 500 Myrs. The star 
formation rate is assumed to be proportional to the mass 
of gas to the power 1.5, with a star formation efficiency v s 
of 3, 1 and 0.5 Gyr -1 , respectively, according to Table 1 in 
Mouhcine & Contini (2002) and their equation (3), ip(t) = 
c s M^ for the starbursting phase. The time for continuous 
infall of primordial gas, ti„f a u, is assumed to be 4 Gyrs. 
Figure 9 shows, similar to Figure 8, the metallicity pre- 
dicted by Pegase2 models with these starbursting SFRs for 
a look backtime of 5 Gyrs and for today, compared to the 
observed metallicity-luminosity relations. Different from 
Figure 8 we show here the metallicity-luminosity depen- 
dence found putting together all the emission line galaxies 
at z <~ 0.4 and z ~ 0.64 (dashed line in Figure 9 labeled 
as 5 Gyrs ago). This line corresponds to a metallicity- 
luminosity relation: 

12 + log(0/H) = (3.55 ± 0.60) - (0.251 ± 0.031)M B , (10) 

with an RMS scatter of 0.227. 

Comparison between Figures 8 and 9 shows that the 
metallicity-luminosity relation at medium redshift is dis- 
placed to lower abundances and higher luminosities com- 
pared to today. This can be explained independent of 
the assumption of a starbursting star formation history 
in the Pegase2 models or of an exponentially decreas- 
ing star formation rate. The shift we see in the Figures 
8 and 9 is comparable to the internal scatter of the 



metallicity-luminosity relation today and at medium red- 
shift. However, our emission line sample is not biased com- 
pared to the Melbourne and Salzer local sample. Therefore 
we think that the shift in our mean relation reflects a true 
shift, in the sense that the metallicity of galaxies increases 
and their luminosity decreases between z ~ 0.7 and z <~ 0. 

5.2. Discussion 

Galaxies over a large range of absolute magnitude (—17 ^ 
Mb Si —22) show an evolution of the metallicity- 
luminosity at z ~ 0.4 and z ~ 0.64 compared to z = 0. 
Comparison with Pcgase2 models shows a possible sce- 
nario in which galaxies at medium redshift have faded by 
~ 0.5 mag (galaxies with lower luminosities in Fig. 8) up 
to <~ 0.9 mag (galaxies with higher luminosities in Fig. 8) 
in <~ 6 Gyrs due to decreasing levels of star formation, and 
their metallicity has increased by a factor of ~ 0.3 dex. 

The luminosity evolution of AM B = 0.5 — 0.9 mag be- 
tween z = and z = 0.64 (5.9 Gyr ago) is consistent 
with the observational results of the COMBO-17 survey 
(Wolf et al. 2003, their Fig.17): They find that M B of star- 
burst galaxies decreases by ~ 0.5 mag between z = 0.7 and 
z = 0.3. 

Fig. 8 and 9 show an increase of metallicity by a fac- 
tor of ~ 2 for galaxies in the local universe compared to 
medium redshift. The metallicity evolution is consistent 
with models of Somerville & Primack (1999) and Pei et al. 
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Fig. 9. The metallicity-luminosity relation, similar to Figure 8, but for starbursts models of Mouhcine & Contini (2002). 
Moreover, the dashed line labeled as 5 Gyrs look back time is the metallicity-luminosity relation found puting together the 
metallicities of all emission line galaxies at z ~ 0.4 and z ~ 0.64. 



(1999), who predict a change in the star-forming metallic- 
ity with redshift of 0.2 dex, and 0.3 dex, respectively, over 
the last half of the age of the universe (see also Fig. 18 
in Lilly et al. 2003). Our result is also consistent with the 
finding of Lilly et al. (2003): Based on estimates of the 
oxygen abundance in a sample of 66 CFRS galaxies at 
0.47 < z < 0.92 they argue that luminous low metallicity 
galaxies at medium redshift, which have normal mass-to- 
light ratios and relatively large sizes, will evolve to higher 
metallicity galaxies of similar absolute magnitudes rather 
than to fainter low metallicity galaxies. 
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